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Abstract 

We consider the Schrodinger operator on graphs and study the spectral statis- 
tics of a unitary operator which represents the quantum evolution, or a quan- 
tum map on the graph. This operator is the quantum analogue of the classical 
evolution operator of the corresponding classical dynamics on the same graph. 
We derive a trace formula, which expresses the spectral density of the quan- 
tum operator in terms of periodic orbits on the graph, and show that one 
can reduce the computation of the two-point spectral correlation function to 
a well defined combinatorial problem. We illustrate this approach by con- 
sidering an ensemble of simple graphs. We prove by a direct computation 
that the two-point correlation function coincides with the CUE expression for 
2x2 matrices. We derive the same result using the periodic orbit approach 
in its combinatorial guise. This involves the use of advanced combinatorial 
techniques which we explain. 
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I. INTRODUCTION 



We have recently shown [1110 that the Schrddinger operator on graphs provides a useful 
paradigm for the study of spectral statistics and their relations to periodic orbit theory. In 
particular, the universal features which are observed in quantum systems whose classical 
counterpart is chaotic, appear also in the spectra of quantum graphs. This observation 
was substantiated by several numerical studies. The relevance to quantum chaology was 
established by identifying the underlying mixing classical evolution on the graphs, which 
provides the stability coefficients and actions of periodic orbits in whose terms an exact 
trace formula can be written ||[IJ,|2|] . 

In spite of the large amount of effort invested in the past fifteen years 0||| , we have only 
a limited understanding of the reasons for the universality of spectral statistics in systems 
whose classical dynamics is chaotic. The main stumbling block is the lack of understanding 
of the intricate and delicate interference between the contributions of (exponentially many) 
periodic orbits. This genuinely quantum quantity, (also known as the "off-diagonal" con- 
tribution), is the subject of several researches, which address it from various points of view 
. The present contribution attempts to illuminate this issue from yet another angle, 
and we harness for this purpose quantum graphs and combinatorics. 

Our material is presented in the following way. We shall start by defining the quantum 
dynamics on the graph in terms of a quantum map. This map will be represented by a 
unitary matrix, which is the quantum analogue of the classical Frobenius-Perron operator 
of the properly defined classical dynamics on the graph. The spectrum of the quantum 
operator is on the unit circle, and its statistics is the main object of the present work. After 
defining the two-point correlation function of interest, we shall write it down in terms of 
periodic orbits and discuss the combinatorial problem which should be addressed in order to 
obtain a complete expression which includes the "off-diagonal" contribution. Since the RMT 
is known to reproduce the two-point correlation function for generic graphs, we propose that 
the RMT expression could be obtained from a combinatorial theory, perhaps as the leading 
term in an asymptotic expansion. For one particular example we show that this is indeed the 
case in the last section. There we construct an ensemble of simple graphs with non-trivial 
spectral statistics, which can be solved in two independent ways. The direct way yields the 
statistics of RMT for the 2x2 circular unitary ensemble (CUE). The corresponding periodic 
orbit calculation is converted into a combinatorial problem, which is solved by proving a 
previously unknown combinatorial identity. 



II. THE QUANTUM SCATTERING MAP AND ITS CLASSICAL ANALOGUE 

A. General Definitions for Quantum Graphs 

We shall start with a few general definitions. Graphs consist of V vertices connected by 
B bonds (or edges). The valency v% of a vertex % is the number of bonds meeting at that 
vertex. Associated to every graph is its connectivity (adjacency) matrix Cjj. It is a square 
matrix of size V whose matrix elements Cij are given in the following way 

„ „ [1 if i, ? are connected I ,. . T . N /1N 

r "' ( - {(. otherwise ) (hJ = l,-,V). (I) 
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The valency of a vertex is given in terms of the connectivity matrix, by t>j = J2j=i an d 
the total number of bonds is B = | Sij=i Cij- 

When the vertices i and j are connected, we shall assume that the connection is achieved 
by a single bond, such that multiple bonds are excluded. We denote the connecting bond by 
b = Note that the notation will be used whenever we do not need to specify the 

direction on the bond. Hence = Directed bonds will be denoted by and we 

shall always use the convention that the bond is directed from the first index to the second 
one. To each bond we assign a length = = Lq^. In most applications we 
would avoid non-generic degeneracies by assuming that the Lfy] are rationally independent. 
The mean length is defined by (L) = -g J2b=i Lb- 

For the quantum description we assign to each bond b = a coordinate Xb which 
measures distances along the bond. We may use x^j) which is defined to take the value at 
the vertex i and the value Luj\ = at the vertex j. We can also use xun which vanishes 
at j and takes the value Lu t j\ at i. 

The wave function $ is a B— component vector and will be written as (^>b 1 {xb 1 ), 
^62(^62)) • • •> ^b B (xb B )) T where the set {h}f =1 consists of all the B distinct bonds on the 
graph. We will call *$?b(xb) the component of \I/ on the bond b. The bond coordinates Xb were 
defined above. When there is no danger of confusion, we shall use the shorthand notation 
*$?b(x) for ^b{xb) and it is understood that x is the coordinate on the bond b to which the 
component ^b refers. 



The Schrodinger equation is defined on the graph in the following way [|H)i|iT| (see also 
|| for an extensive list of references on the subject): On each bond b, the component ^ of 
the total wave function \I/ is a solution of the one-dimensional equation 



(-i d/dx^j) - ty b (x(i,j)) = k 2 ^ b (x(i,j)) (b = ■ (2) 

We included a "magnetic vector potential" Auj\, with Auj\ = —Au^ which breaks time- 
reversal symmetry. 

On each of the bonds, the general solution of (0) is a superposition of two counter- 
propagating waves 

^(i,j)( X (i,j)) = eX P kx (i,j) + A {i,j) X {i,j) ) 

^U,i)( X U,i)) = eX P kx U,i) + A U,i) X U,i) ) • ( 3 ) 

Note that the above functions are normalised to have an amplitude 1 at the points from 
which they "emerge", namely, ip(%j) = 1 at the vertex i and V'Om) = 1 at the vertex j. The 
Hilbert space of the solutions of (0) is spanned by the set of functions defined above, such 
that for all b = 

^b = a(i,j)i>(i,j)(x {iJ) ) + a Uti) tfj {jti) {x {jti) ) . (4) 

Thus, the yet undetermined coefficients au j) form a 25-dimensional vector of complex num- 
bers, which uniquely determines an element in the Hilbert space of solutions. This space 
corresponds to "free wave" solutions since we did not yet impose any conditions which the 
solutions of (H) have to satisfy at the vertices. 
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B. The Quantum Scattering Map 



The quantum scattering map is a unitary transformation acting in the space of free waves, 
and it is defined as follows. 

In a first step, we prescribe at each vertex i — 1, . . . , V a vertex scattering matrix which 
is a unitary matrix of dimension v j. The vertex scattering matrices may be k dependent and 
they are denoted by cr^(/c), where the indices I, m take the values of the vertices which are 
connected to i, that is, C it i = Cj )Tn = 1. The vertex scattering matrix is a property which 
is attributed to the vertex under consideration. It can either be derived from appropriate 
boundary conditions as in |I]|2], or, it can be constructed to model other physical situations. 
The important property of oy^fc) in the present context is, that any wave which is incoming 
to the vertex i from the bonds (I, i), and which has an amplitude 1 at the vertex, is scattered 
and forms outgoing waves in the bonds (i,m) with amplitudes cr^(fc). 

Now, the quantum scattering map is represented by its effect on the 25-dimensional 
vector of coefficients a = ja^j, namely, a is mapped to a' with components 

IB 

a 'b' = X/ a b S B b b , , (5) 

6=1 

where b and b' run over all directed bonds, and if we denote b = and b' = (I, m) 

SB {i , m , m) {k) = Sjj^^^a^k) • (6) 

The effect of Sb on a wave function can be intuitively understood as follows. The coefficient 
au t j\ is the (complex) amplitude of the wave which emerges from the vertex i and propagates 
to the vertex j. Once it reaches the vertex j, it has accumulated a phase e* L ( iJ ^ fc+A ( iJ ^ and 
it scatters into the bonds which emanate from j with an amplitude given by the appropriate 
vertex scattering matrix. The new amplitude a/j =J - m \ consists of the superposition of all the 
amplitudes contributed by waves which impinge on the vertex / = j and then scatter. The 
name "quantum scattering" map is justified by this intuitive picture. 

The resulting matrix Sb is a 2B x IB unitary matrix. The unitarity follows simply from 
the unitarity of the vertex scattering matrices, and from the fact that Sb has non-vanishing 
entries between connected directed bonds: the incoming bond aims at the vertex from which 
the outgoing bond emerges. The unitarity of Sb implies that its spectrum is restricted to 
the unit circle. In this paper we shall mainly be concerned with the spectral statistics of 
the eigenphases, and their relation to the underlying classical dynamics on the graph. The 
spectral statistics will be discussed in the next chapter. We shall use the remaining part of 
the present chapter to clarify two important issues. We shall first show how one can use the 
quantum scattering map to construct the space of solutions of the Schrodinger operator on 
the graph with boundary conditions. Then, we shall introduce the classical dynamics which 
corresponds to the scattering map. 

To define the space of "bound states" on the graph, one has to restrict the space of 
wave functions by imposing appropriate boundary conditions on the vertices. The bound- 
ary conditions guarantee that the resulting Schrodinger operator is self-adjoint. In ]IJ,3, 
we described and used one particular set of boundary conditions, which ensure continu- 
ity (uniqueness) and current conservation. Here we shall use a slight generalisation, which 
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matches well with the spirit of the present article. We shall impose the boundary conditions 
in terms of a consistency requirement that the coefficients auj) have to obey. Namely, we 
require that the wave function (|J) is stationary under the action of the quantum scattering 
map. In other words, the vector a must be an eigenvector of Ssik) with a unit eigenvalue, 
(see also ||12|| ). This requirement can be fulfilled when 

det{I - S B {k)) =0. (7) 

In |Ty2j we have actually derived (0), for the particular case in which the vertex scattering 
matrices where computed form a particular set of vertex boundary conditions which impose 
continuity and current conservation on the vertices. The resulting vertex scattering matrices 
read 

((1 ~\~ e ~ ' lu>i ) \ \ 
—Sjj' H I CijCij', uoi = 2 arctan — l - . (8) 
Vi ) Vik 

Here, < \ < oo are arbitrary constants. The "Dirichlet" ("Neumann") boundary condi- 
tions correspond to Aj = oo (0), respectively. The Dirichlet case implies total reflection at 
the vertex, aj 1 ], = —b~jj>- For the Neumann boundary condition we have <rj l j, = —Sjji + 2/vi 
which is independent of k. For any intermediate boundary condition, the scattering matrix 
approaches the Neumann expression as k — > oo. Note that in all non-trivial cases (v^ > 2), 
back-scattering (j = j') is singled out both in sign and in magnitude: afj has always a neg- 
ative real part, and the reflection probability |o"jj| 2 approaches 1 as the valency i>j increases. 
One can easily check that is a symmetric unitary matrix, ensuring flux conservation 
and time reversal symmetry at the vertex. For Neumann boundary conditions is a real 
orthogonal matrix. 

The spectral theory of the Schrodinger operators on graphs can be developed using (|7J) 
as the starting point. In particular, the corresponding trace formula || can naturally be 
derived, and related to the underlying classical dynamics Here, we shall study the 

quantum scattering map on its own right, without a particular reference to its role in the 
construction of the spectrum. We shall consider the ensemble of unitary, 2B x 2B matrices 
Ssik), where k is allowed to vary in a certain interval to be specified later. Our main concern 
will be the statistical properties of the eigenvalues of Sb- This will be explained in the next 
chapter. 



C. The Classical Scattering Map 

The last point to be introduced and discussed in the present chapter is the classical 
dynamics on the graph and the corresponding scattering map. 

We consider a classical particle which moves freely as long as it is on a bond. The 
vertices are singular points, and it is not possible to write down the analogue of Newton's 
equations at the vertices. Instead, one can employ a Liouvillian approach based on the study 
of the evolution of phase-space densities. This phase-space description will be constructed 
on a Poincare section which is defined in the following way. Crossing of the section is 
registered as the particle encounters a vertex, thus the "coordinate" on the section is the 
vertex label. The corresponding "momentum" is the direction in which the particle moves 
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when it emerges from the vertex. This is completely specified by the label of the next vertex 
to be encountered. In other words, 



position 1 J vertex index 

momentum f | next index 



(9) 



The set of all possible vertices and directions is equivalent to the set of 2B directed bonds. 
The evolution on this Poincare section is well defined once we postulate the transition 
probabilities PfXy between the directed bonds 6 = {j,i} and b' = {i,f}- To make the 
connection with the quantum description, we adopt the quantum transition probabilities, 
expressed as the absolute squares of the Sb matrix elements 



p(i) 



3^3 



" • (io) 



When the vertex scattering matrices are constructed from the standard matching conditions 
on the vertices (H), we get the explicit expression 



p (0 



3~*3 



(1 + e-"* 



(11) 



For the two extreme cases corresponding to Neumann and Dirichlet boundary conditions 
this results in 

(i) f (Sjf + 2/vif Neumann ) , . 

W \ 5 jtj < Dirichlet J ' { ' 

The transition probability Pj%j> for the Dirichlet case admits the following physical inter- 
pretation. The particle is confined to the bond where it started and thus the phase space 
is divided into non-overlapping ergodic components (~ "tori"). For all other boundary 
conditions the graph is dynamically connected. 

The classical Frobenius-Perron evolution operator is a 2B x 2B matrix whose elements 
Ub,y are the classical transition probabilities between the bonds b, b' 

jr. . -a. M'Vi 

U does not involve any metric information on the graph, and for Dirichlet or Neumann 
boundary conditions U is independent of k. This operator is the classical analogue of the 
quantum scattering matrix Sb- Usually, one "quantises" the classical operator to generate 
the quantum analogue. For graphs the process is reversed, and the classical evolution is 
derived from the more fundamental quantum dynamics. 

Let Pb(t), 6=1,..., 2B denote the distribution of probabilities to occupy the directed 
bonds at the (topological) time t. This distribution will evolve after the first return to the 
Poincare section according to 

Pb {t + l)=J2 U b,vPb'{t). (14) 
b' 

This is a Markovian master equation which governs the evolution of the classical probability 
distribution. The unitarity of the graph scattering matrix Sb guarantees J2l=i Ub,v — 1 an d 
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< Ub,b' < 1, such that the probability that the particle is on any of the bonds is conserved 
during the evolution. The spectrum of U is restricted to the unit circle and its interior, and 
vi — 1 is always an eigenvalue with the corresponding eigenvector |1) = ^ (1, 1, 1) T . In 
most cases, the eigenvalue 1 is the only eigenvalue on the unit circle. Then, the evolution 
is ergodic since any initial density will evolve to the eigenvector 1 1) which corresponds to a 
uniform distribution (equilibrium). 

p(t)t=^\l). (15) 

The mixing rate — In | z^2 1 at which equilibrium is approached is determined by the gap 
between the next largest eigenvalue u 2 and 1. This is characteristic of a classically mixing 
system. 

However, there are some non-generic cases such as, e.g., bipartite graphs when —1 belongs 
to the spectrum. In this case the asymptotic distribution is not stationary. Nevertheless an 
equivalent description is possible for bipartite graphs when U is replaced by U 2 which has 
then two uncoupled blocks of dimension B. The example that we are going to discuss in the 
last section will be of this type. 

Periodic orbits on the graph will play an important role in the sequel and we define them 
in the following way. An orbit on the graph is an itinerary (finite or infinite) of successively 
connected directed bonds {^1,^2}, {^2, h}, • • • For graphs without loops or multiple bonds this 
is uniquely defined by the sequence of vertices 11,12, ■ ■ ■ with i m e [1, V] and Cj mi j m+1 = 1 for 
all m. An orbit is periodic with period n if for all k, (i n+ k, i n +k+i) = («fc, ik+i)- The code of a 
periodic orbit of period n is the sequence of n vertices i±, . . . , i n and the orbit consists of the 
bonds (im,i m +i) (with the identification i m+n = i m ). In this way, any cyclic permutation of 
the code defines the same periodic orbit. 

The periodic orbits (PO's) can be classified in the following way: 

• Irreducible periodic orbits - PO's which do not intersect themselves such that any 
vertex label in the code can appear at most once. Since the graphs are finite, the 
maximum period of irreducible PO's is V. To each irreducible PO corresponds its 
time reversed partner whose code is read in the reverse order. The only PO's which 
are both irreducible and conjugate to itself under time reversal are the PO's of period 
2. 

• Reducible periodic orbits - PO's whose code is constructed by inserting the code of any 
number of irreducible PO's at any position which is consistent with the connectivity 
matrix. All the PO's of period n > V are reducible. 

• Primitive periodic orbits - PO's whose code cannot be written down as a repetition of 
a shorter code. 

We introduced above the concept of orbits on the graph as strings of vertex labels whose 
ordering obeys the required connectivity. This is a finite coding which is governed by a 
Markovian grammar provided by the connectivity matrix. In this sense, the symbolic dy- 
namics on the graph is Bernoulli. This property adds another piece of evidence to the 
assertion that the dynamics on the graph is chaotic. In particular, one can obtain the 
topological entropy Y from the symbolic code. Using the relation 
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r = lim -logtr(C n ) (16) 

one gets T = logv, where v is the mean valency. 

Of prime importance in the discussion of the relation between the classical and the 
quantum dynamics are the traces u n = tr(U n ) which are interpreted as the mean classi- 
cal probability to perform n-periodic motion. Using the definition ([13]) one can write the 
expression for u n as a sum over contributions of n-periodic orbits 

u n = n p exp(-r7 p n p ) , (17) 

where the sum is over the set V n of primitive PO's whose period n p is a divisor of n, 
with r = n/n p . To each primitive orbit one can assign a stability factor exp(— 7 P n p ) which 
is accumulated as a product of the transition probabilities as the trajectory traverses its 
successive vertices: 

rip 

eM-%n P ) = II • ( 18 ) 

The stability exponents 7 P correspond to the Lyapunov exponents in periodic orbit theory. 

When only one eigenvalue of the classical evolution operator U is on the unit circle, one 
has, u n n=^ 1. This leads to a classical sum-rule 

u n = n p exp(— r"fpn p ) n=^ 1. (19) 
peP„ 

This last relation shows again that the number of periodic orbits must increase exponentially 
with n to balance the exponentially decreasing stability factors of the individual periodic 
orbits. The topological entropy can be related to the mean stability exponent through this 
relation. 

Using the expression (|T7|) for u n one can easily write down the complete thermodynamic 
formalism for the graph. Here, we shall only quote the periodic orbit expression for the 
Ruelle ( function 



C R (z) = (det(J - zU))~ l = exp [-tr (ln(J - zU))} (20) 

.TT I , 

V (1 - z n v exp(-n P 7 p )) ' 



exp 



„ n 



where the product extends over all primitive periodic orbits. 

The above discussion of the classical dynamics on the graph shows that it bears a striking 
similarity to the dynamics induced by area preserving hyperbolic maps. The reason under- 
lying this similarity is that even though the graph is a genuinely one-dimensional system, 
it is not simply connected, and the complex connectivity is the origin and reason for the 
classically chaotic dynamics. 
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III. THE SPECTRAL STATISTICS OF THE QUANTUM SCATTERING MAP 



We consider the matrices Sb defined in (^). Their spectrum consist of IB points con- 
fined to the unit circle (eigenphases) . Unitary matrices of this type are frequently studied 
since they are the quantum analogues of classical, area preserving maps. Their spectral 
fluctuations depend on the nature of the underlying classical dynamics |13| . The quantum 



analogues of classically integrable maps display Poissonian statistics while in the opposite 
case of classically chaotic maps, the statistics of eigenphases conform quite accurately with 
the results of Dyson's random matrix theory (RMT) for the circular ensembles. The ensem- 
ble of unitary matrices which will be used for the statistical study will be the set of matrices 
Ssik) with k in the range \k — ko\ < A&/2. The interval size A& must be sufficiently small 
such that the vertex matrices do not vary appreciably when k scans this range of values. 
Then the k averaging can be performed with the vertex scattering matrices replaced by their 
value at k . When the vertex scattering matrices are derived from Neumann or Dirichlet 
boundary conditions, the averaging interval is unrestricted because the dimension of Sb is 
independent of k. In any case A^ must be much larger than the correlation length between 
the matrices Ss{k), which was estimated in M to be inversely proportional to the width of 
the distribution of the bond lengths. The ensemble average with respect to k will be denoted 
by 

1 rk +A k /2 

(,) *= wL a „ - dk - (21) 

A fe Jk -A k /2 

Another way to generate an ensemble of matrices Sb is to randomise the length matrix L or 
the magnetic vector potentials ^.(ij), while the connectivity (topology of the graph) is kept 
constant. In most cases, the ensembles generated in this way will be equivalent. In the last 
section we will also consider an additional average over the vertex scattering matrices. 

In the following subsections we compare statistical properties of the eigenphases {Oi(k)} 
of Sb with the predictions of RMT |~7| and with the results of periodic orbit theory for 
the spectral fluctuations of quantised maps [ 18[| . The statistical measure which we shall 



investigate is the spectral form factor. Explicit expressions for this quantity are given by 
RMT |L4| , and a semiclassical discussion can be found in @,|(],|16[ . 



A. The Form Factor 

The matrix Sb for a fixed value of k is a unitary matrix with eigenvalues e l9l( - k \ The 
spectral density of the eigenphases reads 

2B r>r5 i oo 

d(6; k) = £ 6(9 - 0,(*O) = «- + 7T E ^ ien ^S n B {k) + c.c. , (22) 

1=1 271 271 n=l 

where the first term on the r.h.s. is the smooth density d — |p. The oscillatory part is a 
Fourier series with the coefficients tr5*g(/c). This set of coefficients will play an important 
role in the following. Using the definitions (§) one can expand trS l ^(/c) directly as a sum 
over n— periodic orbits on the graph 
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tiS%(k) = 53 n p A r / {kl ^ )r e illpr , (23) 

where the sum is over the set V n of primitive PO's whose period n p is a divisor of n, with 
r = n/n p . l p = Ylbep Lb is the length of the periodic orbit. $ p = J2b&p ^b-^-b is the "magnetic 
flux" through the orbit. If all the parameters Ab have the same absolute size A we can write 
$ p = Ab p , where b p is the directed length of the orbit. \x v is the phase accumulated from 
the vertex matrix elements along the orbit, and it is the analogue of the Maslov index. For 
the standard vertex matrices (§) /x p /7r gives the number of backs catterings along p. The 
amplitudes A p are given by 

rip 

A = nhfti^+i|^-^, (24) 

where ij runs over the vertex indices of the periodic orbit, and j is understood modrip. 
The Lyapunov exponent 7 P was defined in fll8l) . It should be mentioned that (|23|) is the 
building block of the periodic orbit expression for the spectral density of the graph, which 
can be obtained starting from the secular equation (^). In the quantisation of classical area 
preserving maps similar expressions appear as the leading semiclassical approximations. In 
the present context (p3|) is an identity. 

The two-point correlations are expressed in terms of the excess probability density i?2( r ) 
of finding two phases at a distance r, where r is measured in units of the mean spacing |j 

^■M-it r i^)^Sl\\. (25) 

The form factor 

K(n/2B) = ^< \trS n B \ 2 > k (26) 

is the Fourier transform of R2{r, ho). For a Poisson spectrum, K(n/2B) = 1 for all n. RMT 
predicts that K(n/2B), depends on the scaled time n/2B only ||13fl , and explicit expressions 
for the orthogonal and the unitary circular ensembles are known fl~4}| . 

As was indicated above, if the vertex scattering matrices are chosen by imposing Dirichlet 
boundary conditions on the vertices, the classical dynamics is "integrable" . One expects 
therefore the spectral statistics to be Poissonian, 

K(n/2B) = 1 for all n > 1 . (27) 

For Dirichlet boundary conditions the vertex scattering matrices (§) couple only time re- 
versed bonds. Sb is reduced to a block diagonal form where each bond and its time reversed 
partner are coupled by a 2 x 2 matrix of the form 

S (b \k,A)={%^ A)Lb e ). (28) 



c 



The spectrum of each block is the pair ±e lkLb , with the corresponding symmetric and anti- 
symmetric eigenvectors ±1). As a result, we get 
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K(n/2B) = 1 + (-l) n for all n > 1 



(29) 



This deviation from the expected Poissonian result is due to the fact that the extra sym- 
metry reduces the matrix Sb further into the symmetric and antisymmetric subspaces. The 
spectrum in each of them is Poissonian, but when combined together, the fact that the 
eigenvalues in the two spectra differ only by a sign leads to the anomaly (E3). 
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FIG. 1. Form factor for a fully connected graph with V = 20 (a) with and (b) without 
time-reversal symmetry. The smooth curves show the predictions of the corresponding random 
matrix ensembles COE and CUE, respectively. 

Having successfully disposed of the integrable case, we address now the more general 
situation. In Fig. [I] we show typical examples of form factors, computed numerically for a 
fully connected graph with V = 20. The data for Neumann boundary conditions and A = 
(Fig. 1(a)) or A ^ (Fig. 1(b)) are reproduced quite well by the predictions of RMT, which 
are shown by the smooth lines. For this purpose, one has to scale the topological time n by 
the corresponding "Heisenberg time" which is the dimension of the matrix, i.e., 2B. The 
deviations from the smooth curves are not statistical, and cannot be ironed out by further 
averaging. Rather, they are due to the fact that the graph is a dynamical system which 
cannot be described by RMT in all detail. To study this point in depth we shall express the 
form factor in terms of the PO expression (E3). 
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2 



K(n/2B) 



1 



r i(klp+Ab p +iTfi p )r 



(30) 



2B 



k 



1 



n p n p/ A r p Ap, exp {iA(rb p - r'6 p/ ) + in(rfi p - r'/i P /)} 
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The k averaging is carried out on such a large interval that the double sum above is restricted 
to pairs of periodic orbits which have exactly the same length. The fact that we choose the 
lengths of the bonds to be rationally independent will enter the considerations which follow 
in a crucial way. 

The largest deviations between the numerical data and the predictions of RMT occur for 
n — 1, 2. For n = 1 one gets instead of the COE (CUE) values 1/B (1/2B), simply because 
the graph has no periodic orbits of period 1. This could be modified by allowing loops, which 
were excluded here from the outset. The 2-periodic orbits are self-retracing (i.e. invariant 
under time reversal), and each has a distinct length. Their contribution is enhanced because 
back scattering is favoured when the valency is large. Self-retracing implies also that their 
contribution is insensitive to the value of A. The form factor for n = 2 calculated for a fully 
connected graph with v = V — 1 is 



independent of the value of A. This is different from the value expected from RMT. The 
repetitions of the 2-periodic orbits are also the reason for the odd-even staggering which is 
seen for low values of r = n/2B. They contribute a term which is ~ 2exp(— 2Vt) and thus 
decays faster with the scaled time r when the graph increases. 

The deviations between the predictions of RMT and periodic orbit theory for low values 
of r are typical and express the fact that for deterministic systems in general, the short time 
dynamics is not fully chaotic. The short time domain becomes less prominent as B becomes 
larger because the time n has to be scaled by 2B. This limit is the analogue of the limit 
h — > in a general system. 

Consider now the domain 2 < n <C 2B. The PO's are mostly of the irreducible type, 
and the length restriction limits the sum to pairs of orbits which are conjugate under time 
reversal. Neglecting the contributions from repetitions and from self-retracing orbits we get 



The classical return probability u n approaches 1 as n increases (see (|19|) ). Neglecting the 
short time deviations, we can replace u n by 1, and we see that the remaining expression is the 
classical expectation of cos 2 Ab p over PO's of length n. For A = this factor is identically 
1 and one obtains the leading term of the COE expression for n 2B. If A is sufficiently 
large (cos 2 Ab p ) n ~ 1/2, one obtains the short-time limit of the CUE result. The transition 
between the two extreme situations is well described by 




(31) 




(32) 




(33) 
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This formula is derived by assuming that the total directed length b p of a periodic orbit is 
a sum of elementary lengths with random signs. 

The basic approximation so far was to neglect the interference between contributions of 
periodic orbits with different codes (up to time reversal). This can be justified as long as 
periodic orbits with different codes have different lengths. This is the case for low values of 
n. As n approaches B the degeneracy of the length spectrum increases, and for n > 2B all 
the orbits are degenerate. In other words, the restriction rl p = r'l p i in (RO) does not pick up 



a unique orbit and its time reversed partner, but rather a group of isometric but distinct 
orbits. Therefore, the interference of the contributions from these orbits must be calculated. 
The relative sign of the terms is determined by the "Maslov" index. The computation 
of the interfering contributions from different periodic orbits with neighbouring actions is 
an endemic problem in the semiclassical theory of spectral statistics. These contributions 
are referred to as the non-diagonal terms, and they are treated by invoking the concept of 
periodic orbit correlations The dynamical origin of these correlations is not known. 

In the case of graphs, they appear as correlations of the "Maslov" signs within a class of 
isometric n-periodic orbits. 

To compute K{n/2B) from fl3~0|) one has to sum the contributions of all the n-periodic 
orbits after grouping together those which have exactly the same lengths. We shall discuss 
the case A = 0, so a further restriction on the orbits to have the same directed length is 
not required here. Since the lengths of the individual bonds are assumed to be rationally 
independent, a group of isometric n-periodic orbits is identified by the non-negative integers 
qi, i = 1, . . . , B such that 

B B 

l 1 = ^2 Qih witn E Qi = n > ( 34 ) 
1=1 1=1 

i.e., each bond i is traversed times. The orbits in the group differ only in the order by which 
the bonds are traversed. We shall denote the number of isometric periodic orbits by D n (q). 



Note that not all the integer vectors q which satisfy (34) correspond to periodic orbits. 



Rather, the connectivity required by the concept of an orbit imposes restrictions, which 



render the problem of computing D n (q) a very hard combinatorial problem [15|. Writing 
(j30l) explicitly for the case of a fully connected graph with Neumann vertex scattering 
matrices, we get 



K(n/2B) = — (-) £ 



A.(q) .. 

E -(-0' 



with f = — — , (35) 



and the a summation extends over the n-periodic orbits in the class q. u Q is the number of 
back scattering along the orbit, and r a is different from unity if the orbit is a repetition of 
a shorter primitive orbit of period n/r a . 



Equation ( |35"D is the starting point of the new approach to spectral statistics, which we 
would like to develop in the present paper. The actual computation of (|35|) can be consid- 
ered as a combinatorial problem, since it involves counting of loops on a graph, and adding 
them with appropriate (signed) weights. For Neumann boundary conditions, the weights 
are entirely determined by the connectivity of the graph. Our numerical data convincingly 
show that in the limit of large B the form factors for sufficiently connected graphs repro- 
duce the results of RMT. The question is, if this relation can be derived using asymptotic 
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combinatorial theory. The answer is not yet known, but we would like to show in the next 
section that for a very simple graph one can use combinatorics to evaluate the periodic orbit 
sums, and recover in this way the exact values of the form factor. 



IV. THE 2-STAR MODEL 



In this section we will investigate the classical and quantum dynamics in a very sim- 
ple graph using two different methods. We shall use periodic orbit theory to reduce the 
computation of the trace of the classical evolution operator u n and the spectral form factor 
K(n/2B) to combinatorial problems, namely sums over products of binomial coefficients. 
The result will be compared to a straight forward computation starting from the eigenvalues 
of the classical and quantum scattering maps. 

An n-star graph consists of a "central" vertex (with vertex index 6) out of which emerge 
n bonds, all terminating at vertices (with indices j = 1, . . . , n) with valencies Vj = 1. The 
bond lengths are L j = Lj. This simple model (sometimes called a hydra) was studied 
at some length in ||. The star with n = 2 is not completely trivial if the central vertex 
scattering matrix is chosen as 



where the value 0<?7<7r/2is still to be fixed. The scattering matrices at the two other 
vertices are taken to be 1 and correspond to Neumann boundary conditions. The dimension 
of U and Sb is 4, but it can be immediately reduced to 2: due to the trivial scattering at 
the reflecting tips, for j = 1,2. In this representation the space is labelled 

by the indices of the two loops (of lengths 2L\ and 2L 2 respectively) which start and end at 
the central vertex. After this simplification the matrix Sb reads 



We shall compute the form-factor for two ensembles. The first is defined by a fixed value of 
i] = tt/4, and the average is over an infinitely large k range. The second ensemble includes an 
additional averaging over the parameter t]. We will show that the measure for the integration 
over r) can be chosen such that the model yields the CUE form factor. This is surprising at 
first sight, since the model defined above is clearly time- reversal invariant. However, if we 
replace kL\ and kL 2 in ([37]) by L(k±A), fl3~7|) will allow for an interpretation as the quantum 
scattering map of a graph with a single loop of length L and a vector potential A, i.e., of 
a system with broken time-reversal invariance (see Fig. |2[). In particular, the form factors 
of the two systems will coincide exactly, when an ensemble average over L is performed. 
Clearly, this is a very special feature of the model considered, and we will not discuss it here 
in more detail. 




(36) 




(37) 



A. Periodic Orbit Representation of u. 



"it. 



The classical evolution operator corresponding to ([$7]) is 
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/cos^sinM 
V sin rj cos r\ J 

The spectrum of U consists of {1, cos 27/}, such that 

Un {rj) = 1 + cos n 2r/. (39) 

We will now show how this result can be obtained from a sum over the periodic orbits of the 
system, grouped into classes of isometric orbits. This grouping is not really necessary for a 
classical calculation, but we would like to stress the analogy to the quantum case considered 
below. 

The periodic orbits are uniquely encoded by the loop indices, such that each n-tuple of 
two symbols 1 and 2 corresponds (up to a cyclic permutation) to a single periodic orbit. 
When n is prime, the number of different periodic orbits is N 2 (n) = 2 + (2 n — 2)/n, otherwise 
there are small corrections due to the repetitions of shorter orbits. These corrections are 
the reason why it is more convenient to represent a sum over periodic orbits of length n as 
a sum over all possible code words, though some of these code words are related by a cyclic 
permutation and consequently denote the same orbit. If we do so and moreover replace the 
stability factor of each orbit by ([18]), the periodic orbit expansion of the classical return 
probability becomes 

n 

Un= E £ II , (40) 

h=l,2 i n =l,2j=l 

where j is a cyclic variable such that i n+ \ = i\. In fact fl4"0| ) can be obtained without any 
reference to periodic orbits if one expands the intermediate matrix products contained in 
u n = tiU n and uses Pi s -M j+l = U ijiij+1 (r}). 

We will now order the terms in the multiple sum above according to the classes of 
isometric orbits. In the present case a class is completely specified by the integer q = qi 
which counts the traversals of the loop 1, i.e., the number of symbols 1 in the code word. 
Each of the q symbols 1 in the code is followed by an uninterrupted sequence of tj > 
symbols 2 with the restriction that the total number of symbols 2 is given by 

Y.t^n-q. (41) 

We conclude that each code word in a class < q < n which starts with a symbol i\ — 1 
corresponds to an ordered partition of the number n — q into q non-negative integers, while 
the words starting with i\ = 2 can be viewed as partition of q into n — q summands. 

To make this step very clear, consider the following example: All code words of length 
n = 5 in the class q = 2 are 11222, 12122, 12212, 12221 and 22211, 22121, 21221, 22112, 
21212, 21122. The first four words correspond to the partitions + 3 = 1 + 2 = 2 + 1= 3 + 
of n — q = 3 into q = 2 terms, while the remaining 5 words correspond to 2 = + + 2 = 
+ 1 + 1 = 1 + + 1 = + 2 + = 1 + 1 + = 2 + + 0. 

In the multiple products in fl4"0|), a forward scattering along the orbit is expressed by two 
different consecutive symbols ij ^ in the code and leads to a factor sin 2 1], while a back 
scattering contributes a factor cos 2 rj . Since the sum is over periodic orbits, the number of 
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forward scatterings is always even and we denote it with 2v. It is then easy to see that v 
corresponds to the number of positive terms in the partitions introduced above, since each 
such term corresponds to an uninterrupted sequence of symbols 2 enclosed between two 
symbols 1 or vice versa and thus contributes two forward scatterings. For the codes starting 
with a symbol 1 there are rM ways to choose the v positive terms in the sum of q terms, 

and there are y^^j wa Y s to decompose n — q into v positive summands. After similar 
reasoning for the codes starting with the symbol 2 we find for the periodic orbit expansion 
of the classical return probability 



u n (r ] ) = 2cos 2n V +J2J2 



2cos 2 ^ + ££- 





sin 4 ^ cos 2n ~ Au T] 



■ iu 2n—Av 

sin rj cos rj 




• Av 2n—iv 

sin rj cos r] 



cos T) + sin rf) n + (cos rj — sin r/) 7 



(42) 



which is obviously equivalent to (|39D. The summation limits for the variable v are implicit 
since all terms outside vanish due to the properties of the binomial coefficients. In order to 
get to the third line we have used the identity 

n-l 




1 



It can be derived by some straightforward variable substitutions from 




which, in turn, is found in the literature 




(43) 



(44) 



B. Quantum Mechanics: Spacing Distribution and Form Factor 



Starting from ([$?]), and writing the eigenvalues as e lk ( L ^-+ L ^) e ±lX / 2 ^ we get for A, the 
difference between the eigenphases, 



A = 2 arcos [cos rj cos k(Li — L 2 )] . 



(45) 



For fixed 77, the k averaged spacing distribution (which is essentially equivalent to i? 2 (^) f° r 
the considered model) is given by 

1 /•feo+Afc/2 

P{6;rj) = — / dk 5 (6 — 2arcos [cos^ cos A;(Li — L 2 )D 



A fc Jfco-A fc /2 



sin(0/2) 



^ ^Jcos 2 r]- cos 2 (#/2) 



cos(0/2) > I COS Tj 
cos(#/2) < I cos 77 



(46) 



16 



We have assumed that 6 is the smaller of the intervals between the two eigenphases, i.e. 
< 6 < 7T. 

The spacings are excluded from a domain centered about (it), i.e., they show very 
strong level repulsion. The distribution is square-root singular at the limits of the allowed 
domain. 

P(#; rj) can be written as 

1 1 OO / 1 \ 

P(6 ]V ) = ^ + -E cos(nfl) (j (\tvS B ( V r\ 2 ) k - lj , (47) 
and, by a Fourier transformation, we can compute the form factor 

K 2 (n;r ] ) = ±(\trS B (r 1 r\ 2 ) . (48) 
In particular, for 77 = it/ 4 one finds 

(_l)m+n /o m \ 

/C 2 (,,V4) = l + i i ^ i -( 2 ™) (49) 

( 1 \m+n 

2a/ 7m 

Where m = [n/2] and [■] stands for the integer part. The slow convergence of K2(n; vr/4) to 
the asymptotic value 1 is a consequence of the singularity of P(6>; 7r/4). 

We now consider the ensemble for which the parameter rj is distributed with the mea- 
sure dn(rj) = I cosr/sinr/ld??. The only reason for the choice of this measure is that upon 
integrating (f47|) one gets 



P(6) = 2sin 2 (#/2), (51) 
which coincides with the CUE result for 2x2 matrices. A Fourier transformation results in 

^«)={H:::>2- m 

The form factors (|49|) , (|50|) and ( ^2|) are displayed in Fig. |] below. 



C. Periodic Orbit Expansion of the Form Factor 



As pointed out at the end of section Hi A , the k- averaged form factor can be expressed 
as a sum over classes of isometric periodic orbits. The analogue of (|35D for the 2-star is 



K 2 (n; 77) = -J2 



q=0 



^-^ r 

a=l 1 a 



2;/, 



■2i> a 



1 sin r/cos 



(53) 



where the number of forward and backward scatterings along the orbits are 2v a and /i Q = 
n — 2u a , respectively. Again, it is very inconvenient to work with the repetition number r a , 
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and consequently we replace — as in the derivation of ( 4~2"D — the sum over orbits by a sum 
over all code words and use the analogy with the compositions of integer numbers to obtain 



K 2 {n;r)) = cos 2n r/ + — £ 

z q=l 



2 re— 1 



E 



V 



l) v (q - l\(n-q-l S 



V 



1 



sm 2u r] cos" 2v T] 



(54) 



The inner sum over v can be written in terms of Krawtchouk polynomials [pi, HI 

1 2 



as 



^ n— 1 

^2(71577) = cos 2n r/ + - J2 



n 



i - q. 



" — ,(cos 2 n,sin 2 n) 



-P. 



9 



n— l,re— g 



(55) 



and the Krawtchouk polynomials are defined as in |21,22| by 



N,k 1 



x) 



-1/2 fc 
i/=0 



k—v 



X 



f N - x" 



u k ~ v v v 



< k < N 

U + V = 1 



(56) 



These functions form a complete system of orthogonal polynomials of integer x with < x < 
iV. They have quite diverse applications ranging from the theory of covering codes [23[] to 
the statistical mechanics of polymers [p^j , and are studied extensively in the mathematical 
literature [ 21 , 22 1 . The same functions appear also as a building block in our periodic orbit 
theory of Anderson localisation on graphs |25| . Unfortunately, we were not able to reduce the 
above expression any further by using the known sum-rules and asymptotic representations 
for Krawtchouk polynomials. The main obstacle stems from the fact that in our case the 
three numbers N, k, x in the definition (|56D are constrained byiV = A; + x — 1. 

We will now consider the special case rj = tt/4 for which we obtained in the previous 
subsection the solution (f49|). The result can be expressed in terms of Krawtchouk polynomi- 
als with u = v = 1/2 which is also the most important case for the applications mentioned 
above. We adopt the common practice to omit the superscript (u, v) in this special case and 
find 



1 



It is convenient to introduce 



1 



2 n+l 



n-1 

E 



(n — 


;) 


\n — 





n 



Pn-l,n-q{q) 



(57) 



\s+t 



S + t-1 

s 



1/2 



P. 



s+t-1, 



E(-i)^( 



t 



8-1 



V \V 



(58) 



and to rewrite (57) with the help of some standard transformations of binomial coefficients 

as 

1 2 



1 1 n_1 
K 2 (n;7T/A) = — + -V 

Z Z q=l 



n 



—J\f(q, n — q — 1) 

q 



1 1 



n-1 



J2Mq,n-q) + (-iyAf(n-q,q)Y 



(59) 



9=1 
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This expression is displayed in Fig. || together with fl49|) in order to illustrate the equivalence 
of the two results. An independent proof for this equivalence can be given by comparing the 
generating functions of K 2 {n\ vr/4) in the two representations |27] 



We defer this to appendix 
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FIG. 2. Form factor for the 2-star quantum graph. The crosses and the connecting heavy full 
line show the two equivalent exact results 



and (57) for r] = tt/4. The thin dashed lines represent 
the approximation (^0|), and the thin straight line corresponds to the diagonal approximation, when 
repetitions of primitive periodic orbits are neglected. The heavy dashed line exhibits the form factor 
of a CUE ensemble of 2 x 2 random matrices (|52|), which can be obtained from the 2-star by an 
appropriate averaging over rj. Finally, the inset shows a sketch of the two possible realisations of 
the system: a time-reversal invariant 2-star with bond lengths Li,L 2 or a graph with a single loop 
of length L and a magnetic flux A breaking time-reversal symmetry. 

Please note, that in this way we have found a proof for two identities involving 
Krawtchouk polynomials 

I 2 



2m- 1 




2m 

r; 

Q 



2m— 1,2m— q 



(?) 



}2m+l 




and 



2m 



2m 



^2m + l- g , 



2m + 1 

q 



n 2 



P. 



2m.2m+l 



•)2m+2 



2(-l) 



2m 



m 



(60) 



(61) 



which were obtained by separating even and odd powers of n in (49) and (|57D. To the best 
of our knowledge, ( |60|) and (|BTD were derived here for the first time. 
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Finally we will derive the CUE result (|52|) for the ensemble of graphs defined in the 
previous subsection starting from the periodic orbit expansion (54). We find 

K 2 (n) = r /2 dfi( V )K 2 (n; V ). (62) 
Jo 

Inserting ((54"|), expanding into a double sum and using 

7r/2 d? 7 sin 2 (^ , )+ 1 ? 7 co S 2 ("— / ) +1 r ? = - 1 - ( n ) (63) 
we get 

K 2 (n) = ^ + (64) 
n 2 ^ (-1)"+"' { n \~ 1 fq- l\ (n-q- l\ (q- l\ (n-q- l\ 



— YY 

A(n + l)^^, vv' \v + v' \v-l M ^ - 1 M^'-lA ^ - 1 J 

Comparing this to the equivalent result fl5"2"|) we were again led to a previously unknown 
identity involving a multiple sum over binomial coefficients. It can be expressed as 

S{n,q) = Y,F»y{n,Q) = l i 1 < Q < n) (65) 

with 



In this case, an independent computer-generated proof was found ||29|| , which is based on 
the recursion relation 

q 2 F uy (n, q)-(n-q- l) 2 F u y(n } q + 1) + (n - l)(n - 2g - + 1, g + 1) = . (67) 



This recursion relation was obtained with the help of a Mathematica routine , but it can 
be checked manually in a straight forward calculation. By summing (|6T| ) over the indices 



v, v', the same recursion relation is shown to be valid for S(n, q) |]28|j30[| and the proof is 
completed by demonstrating the validity of ( |65| ) for a few initial values. Having proven ( |65| ) 
we can use it to perform the summation over u, v' in (|64]) and find 



K 2 {n) = -±- + £ = -—— H T _ t(^ — <W) > (68) 

n + 1 ^— J — 1 n + 1 n + 1 

which is now obviously equivalent to the random matrix form factor ( [52] ) . To the best of our 
knowledge, this is the first instance in which a combinatorial approach to random matrix 
theory is employed. 
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V. CONCLUSIONS 



We have shown how within periodic orbit theory the problem of finding the form factor 
(the spectral two-point correlation function) for a quantum graph can be exactly reduced 
to a well-defined combinatorial problem. For this purpose it was necessary to go beyond 
the diagonal approximation and to take into account the correlations between the periodic 
orbits. 

In our model, these correlations are restricted to groups of isometric periodic orbits. 
This fits very well with the results of |7j, where for a completely different system (the Sinai 
billiard), the classical correlations between PO's were analysed and found to be restricted to 
relatively small groups of orbits. The code words of the orbits belonging to one group were 
conjectured to be related by a permutation and a symmetry operation, which is in complete 
analogy to the isometric orbits on graphs. 

Even for the very small and simple graph model that we considered in the last section 
the combinatorial problems involved were highly non-trivial. In fact we encountered pre- 
viously unknown identities which we could not have obtained if it were not for the second 
independent method of computing the form factor. However, since the pioneering work 
documented in []30| the investigation of sums of the type we encountered in this paper is a 
rapidly developing subject, and it can be expected that finding identities like flBUD, ( [61]) and 
( |5"5D will shortly be a matter of computer power. 

The universality of the correlations between periodic orbits in all chaotic systems poses 
the problem to identify the common dynamical reasons for their occurrence and to find 
a common mathematical structure which is capable to describe them. A very interesting 
question in this respect is, if the correlations between PO's in a general chaotic system can 
be related to combinatorial problems. 
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APPENDIX A: PROOF OF EQUIVALENCE FOR EQS. (f^) AND (|^) 

In this appendix we give an independent proof for the equivalence between the two results 
(P)|) and ( fjUp obtained in sections |1V B| and |1V G| , respectively, for the form factor of the 
2-star with 77 = 7r/4. We define the generating function 

00 

G{x) = J2 K M t/4) (2x) n (\x\< 1/2) (Al) 

x=l 
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and find from 



1 1 - 2x 1 1 - Qx 



2 VI + 4a; 2 21 -2x 
On the other hand we have from (1591) 



(A2) 



G(x) = z ^ + G 1 (x) + G 2 (-x) (A3) 

with 

oo 

G 1 (x) = AT 2 (s,t)s s+t (A4) 

s,t=l 

and 

oo 

G 2 (x) = ]T JV(s, *) Af(t, s) x s+t . (A5) 

s,t=l 

A convenient starting point to obtain G\ and G2 is the integral representation 



M(s, t) = / d^ (1 + z-'Y (1 - z)- 1 , (A6) 



where the contour encircles the origin. With the help of ( |A6| ) we find 

g(x,y) = Af(s,t)x s y t 

s,t=l 

1 00 „ 00 



,t=l " s,t=l 

x " 1 1 + Z 



27ri s ^ 7 Z l-a;(l- 



(1 -z) z + y(l + z) 



xy 



(l + y)(l-x + y-2xy) 



(|x|,M<l/V2). (A7) 



The contour |1 + z 1 \ = |1 — z\ = V2 has been chosen such that both geometric series 
converge everywhere on it. Now we have 

1 r A A 1 00 00 

Gl(x2) = (2^ 7 ^ £ ? M{s t t)M{J,H) {xzy{x/zr\ X z')\x/z'f 



X 



s,t=l s',t'=l 

■1 



dzdZ (1 + xz')(l + 4? - ij - 2x 2 zz') (z' + x)(zz' + x[z - z'\ - 2x 2 ) ' ^ 

where \x\ < 1/ \[2 and the contour for z, z' is the unit circle. We perform the double integral 
using the residua inside the contour and obtain 
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Gi(x) = 
In complete analogy we find 

G 2 {x) 

such that 



x 



2x-l 1 VA 



1 — x 



1 Ax 2 + 2x + 1 



2(2x + l) v / 4 



+ 



1 Ax 2 - 2x + 1 



1-z 2x-l VV4x 2 + l l-sy 2 (1 - 2x)VAx 2 + 1 2 



(A9) 



(A10) 



(All) 



The proof is completed by a straightforward verification of the equivalence between the 
rational functions ( |A2|) and ([All ). 
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